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ABSTRACT 

Programs to observe evolution in the M»— a or M,—L relations typically com- 
pare black-hole masses, M», in high-redshift galaxies selected by nuclear activity 
to M. in local galaxies selected by luminosity L, or stellar velocity dispersion a. 
Because AGN luminosity is likely to depend on M., selection effects are different 
for high-redshift and local samples, potentially producing a false signal of evolu- 
tion. This bias arises because cosmic scatter in the M. — a and M. — L relations 
means that the mean log 10 L or log 10 o among galaxies that host a black hole 
of given M,, may be substantially different than the log 10 L or log 10 cr obtained 
from inverting the M, — L or M. — o relations for the same nominal M.. The 
bias is particularly strong at high M., where the luminosity and dispersion func- 
tions of galaxies are falling rapidly. The most massive black holes occur more 
often as rare outliers in galaxies of modest mass than in the even rarer high-mass 
galaxies, which would otherwise be the sole location of such black holes in the 
absence of cosmic scatter. Because of this bias, M, will typically appear to be 
too large in the distant sample for a given L or a. For the largest black holes and 
the largest plausible cosmic scatter, the bias can reach a factor of 3 in M. for 
the M. — a relation and a factor of 9 for the M. — L relation. Unfortunately, the 
actual cosmic scatter is not known well enough to correct for the bias. Measuring 
evolution of the M, and galaxy property relations requires object selection to be 
precisely defined and exactly the same at all redshifts. 
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Subject headings: galaxies: nuclei — Galaxies: Evolution — Galaxies: Funda- 
mental Parameters 



Observing Evolution in the Relationships Between Black Hole Mass and 

Galaxy Properties 

The discovery that most elliptical galaxies and spiral bulges host a black hole at their 



centers, plus the tight relations obs erved between 
ity L or stellar velocity dispersion a ( Dresslerl 



1995 : lMagorrian et al.lll998l : lFerrarese fc Merritt 



1989 



jlack-hole mass M, and galaxy luminos- 



Kormendv 



2000 



1993 



Kormendv fc Ri chstonc 



Gebhardt et al. 2000a; Tremaine et al. 



2002c lHaring fc Rix! 120041 ). suggest that the formation and growth of central black holes is 
deepl y intertwined with that of their host galaxies. Recent theoretical work (e.g.. lHopkins et al 
20061 ) supports this view, and also predicts how the relations between the properties of black 
holes and their host galaxies have changed over time. Direct observation of the evolution 
of the M. — a and M. — L relations over cosmological time would offer unique insight into 
galaxy and black-hole formation. 

There have been many attempts to measure the evolution of the M. relations; we cite 
a partial list of these: 



Shields et al.l (120031 ) examine a sample of quasars at redshifts up to 3.3. The luminosity 



L and velocity dispersion a of the host galaxy in such objects cannot be measured 
reliably because the light from the galaxy is overwhelmed by the quasar flux; instead, 
they use the width of the narrow [O III] emission line as a surrogate for a. They 
estimate the black-hole mass using the "photoionization" method, which is based on 
an empirically calibrated relation involving the continuum luminosity of the active 
galactic nucleus (AGN) and the width of the B./3 or other broad emission lines. The 
theoretical assumptions that underlie this relation are that the bulk velocities of the 
emitting clouds are determined by their orbital motion in the gravitational field of the 
black hole, and that the emitting region is photoionized by an ultraviolet continuum 
spectrum of fixed shape. They find an M. — a relation that is consistent with the local 
one, suggesting that this relation is independent of redshift. 



1 The National Optical Astronomy Observatory is operated by AURA, Inc., under cooperative agreement 
with the National Science Foundation. 
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Using similar meth ods on a larger but low er-redshift quasar sample from the Sloan 
Digital Sky Survey, ISalviander et al.l (120061 ) find that galaxies of a given dispersion at 
z ~ 1 have black- hole masses that are larger by A log 10 M, ~ 0.2 than at z — 0. 



Treu et al.l (120041 ) and lWoo et al.l (120061 ) measured the velocity dispersion of 14 Seyfert 
I galaxies at redshift z ~ 0.36; they used stellar absorption lines, which should provide 
a more direct measure of the dispersion than emission lines. They measured the black- 
hole mass using the photoionization method. They found that galaxies of a given 
dispersion at z ~ 0.36 host black holes having A log 10 M. = 0.62 ± 0.10 ± 0.25. 



Peng et al.l (120061 ) measured the bulge luminosities of 11 quasar hosts in the redshift 
range 1.7 < z < 2.7 using the Hubble Space Telescope (HST), and estimated their 
black-hole masses using the photoionization method. They conclude that the M. — L 
relation at z ~ 2 is close to the relation at z = 0; Alog 10 M. ~ —0.1. This result 
is remarkable since black holes can only grow with time, while elliptical galaxies fade 
with time as their stars die: even if the black holes do not grow at all, passive stellar 
evolution models would predict Alog 10 M. between —0.4 and —0.8 at z ~ 2. If their 
result is correct, then either black holes must be ejected from the galaxy centers and 
replaced with smaller ones, or the galaxy luminosity must grow substantially through 
mergers. 



A common thread among all these investigations and others is that the high-redshift 
sample is selected by some measure of AGN visibility. Black-hole masses and galaxy prop- 
erties are then derived and an M. — o or M. — L relation is fitted to the data. Since the 
low-redshift relations are linear in log 10 M., log 10 L, and log 10 cr (eqs. @] or \5\ below) , the 
high-redshift data are often analyzed by assuming that the slope of the relation is the same 
as at low redshift and estimating the offset in the intercept or zero-point, which may be 
expressed as Alog 10 M. at fixed L or a. Any such offset is then interpreted as evidence that 
the ratio of M. to a or L has changed over time. Of course, the offset can be (and sometimes 
is) equally well expressed as A log 10 L or A log 10 o at fixed black-hole mass. 

Proof of evolution in the M. relations requires demonstrating that the high- and low- 
redshift galaxy samples were assembled with no selection effects that would bias the relation 
between the typical galaxy properties and black- hole masses being measured — or, at least, 
that the selection effects in the high- and low-redshift samples were the same. This may be 
a more difficult task than has commonly been assumed. The local black-hole sample has 
mostly been drawn from "normal" galaxies with quiescent or low levels of nuclear activity. 
Despite the proximity of these galaxies, detecting and weighing their central black holes 
requires exquisite spatial resolution, very high signal-to-noise observations, and elaborate 
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modeling of the observations. Investigating the properties of black holes in nearby galaxies 
remains a frontier problem; to date, dynamically determined black-hole masses are available 
for only slightly more than three dozen galaxies. Since there was little evidence of a black 
hole in most of these galaxies before HST observations were taken (except sometimes for a 
modest rise in the velocity dispersion in spectra taken at ground-based resolution), and since 
black holes are found in almost all nearby galaxies that have been examined carefully with 
HST, the black holes in local quiescent galaxies are selected mainly on the basis of galaxy 
properties such as L or a. 

The techniques used to measure black-hole masses in quiescent galaxies at low redshift 
cannot be used at high redshift, both because the radius of influence of the black hole 
cannot be resolved beyond a few tens of Mpc, and because the galaxies have substantially 
lower surface brightness, by the factor (1 + z) A . Black holes at high redshift are identified 
instead by their association with AGN. The observer thus first locates a black hole that is 
accreting matter and weighs it by the properties of the AGN emission lines (a difficult task, 
but one that we shall not examine in detail here), and then secondly attempts the (still 
difficult) task of measuring the properties of the host galaxy. The existence and properties 
of the AGN depend both on the properties of the black hole (mass, spin, orientation) and on 
the properties of the galaxy (mass inflow rate, orientation, etc.). In short, selection is done 
at low redshift by galaxy properties, and at high redshift by a combination of black-hole and 
galaxy properties that depends on the method used. 

Galaxy samples obtained with different sele ction techni q ues w ill generally satisfy dif- 



ferent M. relations. This effect, analogous to the iMalmquistl (119241 ) bias that is familiar in 
studies of Galactic structure, arises because of the cosmic scatter in the M. relations — there 
is not (so far as we know) an exact 1-1 relation between black-hole mass and any single 
measurable property of galaxies such as L or o. In this paper we argue that determinations 
of the redshift evolution of the M, relations may be strongly biased by selection effects unless 
the same sample selection techniques are used at all redshifts. Correcting for the selection 
bias introduced by the use of different sample selection criteria at different redshifts is ex- 
tremely difficult: to make accurate corrections it is necessary to know both how the selection 
depends on the galaxy and black-hole properties and the cosmic scatter in the M. relations. 
At present we have only crude upper limits to the latter quantity. We shall show that even 
if the cosmic scatter were known, the bias in the samples described above usually can only 
be corrected with additional information, such as how galaxy properties and black-hole mass 
determine the probability distribution of AGN luminosity. 



The selection bias can be especially large at large M., as the following argument shows. 
Consider the M. — L relation. At high L the number density of galaxies falls off rapidly, as 
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illustrated by the steep cutoff in the ISchechterl (119761 ) luminosity function. Cosmic scatter 
in the M. — L relation implies that there is a distribution of black-hole masses at a given L. 
The rare high-mass black holes can arise from either the peak of this distribution in the rare 
galaxies with large L, or from the high-mass tail of the distribution in the more numerous 
galaxies of modest L. If the number density of galaxies is falling off rapidly with L, the 
contribution from galaxies with modest L may actually overwhelm the population of black 
holes of simi lar mass as s ociate d with galaxies of higher L. This problem was first explicitly 
identified by lFine et al.l (120061 ) in the context of t he co rrelation between black- hole mass and 
dark-matter halo mass, and by lSalviander et al. (boo6r i in the context of the M. — a relation, 
but apparently has not been appreciated by most observers studying the evolution of the 
M. relations. 

In this paper we present a more general exposition of the biases incurred when studying 
evolution of the M. relations by comparing samples obtained by different selection criteria 
and at different redshifts. We start by considering the extraction of samples from hypothet- 
ical joint distributions of M. and L or a. We then consider the selection bias that occurs if 
the high-redshift samples are selected by AGN flux or luminosity. We discuss the prospects 
of correcting for selection bias, and argue that selection bias may place fundamental limits 
on the determination of the evolution of the M. relations. We conclude with a brief review of 
attempts to measure the evolution of the M, relations and how they may have been affected 
by selection bias. 



2. The Joint Distribution of Black-Hole Mass and Galaxy Properties 

Understanding object selection bias requires knowledge of the joint probability distri- 
bution of black- hole mass M. and a galaxy property s, which for the present discussion is 
either log 10 L, where L is the rest-frame \^-band galaxy luminosity, or log 10 a, where a is the 
line-of-sight velocity dispersion in the main body of the galaxy (the precise definition of L or 
o does not concern us, so long as it is defined consistently in all samples). The true form of 
the joint distribution unfortunately is unknown; however, we can construct a hypothetical 
form of this distribution that accurately represents the present state of the observations, and 
in any case suffices to show how sample selection bias can occur. 

We write the probability of finding a galaxy in the interval (/i, /i + dfi) and (s, s + ds) 
as u(fi, s) dfids where fi = log 10 (M.). From the definition of conditional probability this can 
be rewritten as 

v(jp,s) = v(p\s)g(s), (1) 
where g{s) ds is the probability that a randomly chosen galaxy in a given volume lies in the 
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interval (s, s + ds) and v(p\s) dp is the probability that the black-hole mass lies in the range 
([/,, p + dp) given that the galaxy property is s. The function g(s) is then given by either the 
volume-limited luminosity function — more properly, the luminosity function of early-type 
galaxy components (ellipticals and spiral bulges), since these are the components that corre- 
late with black-hole mass — or the velocity-dispersion distribution of early-type components. 
Without loss of generality, we may write u(p,s) = h[p — f(s),s] where J h(x,y)dx = 1 
and J xh(x,s)dx = 0, so f(s) is the mean value of p at a given value of the galaxy prop- 
erty s — that is, f(s) is either the M. — L or M. — a relation. The limited observational 
data on black-hole masses in local galaxies is consistent with the hypothesis that h(x, s) is 
independent of s; adopting this hypothesis for simplicity we have 

u(p,s) = h[p-f(s)]g(s). (2) 

The intrinsic variance or "cosmic scatter" in black-hole mass at a given value of galaxy 
property s is a 1 = J x 2 h(x, s)dx. Observational errors may make h(x, s) appear to be yet 
broader, but for this discussion, we assume that observational errors are negligible; the 
selection effects that we are concerned with are solely due to the fact that there is an 
intrinsic and irreducible range of M. at any s. 

It should be stressed that although the functional form (j2j) is consistent with the available 
data, it is far from unique. As a foil, we may consider the form 

v(p,s)=p[s-c(p)]<S>.(p). (3) 

A physical model that motivates f$2§ is one in which the galaxy property s determines the 
black-hole mass p with a cosmic scatter described by the function h, while ([3]) is motivated 
by models in which the black- hole mass p determines the galaxy property s, with cosmic 
scatter described by the function p. If f p(y)dy = 1, then is the number of black 

holes per unit volume with log mass in the range (p, fi + dp). We do not know which of ([2]) 
or ([3]) is correct (possibly neither, or both); our motivation for choosing the former is that it 
provides a simple representation of what local observations of black holes in inactive galaxies 
measure. 

Of the three functions that contribute to u(p,s), q (s) is probab l y the best determined. 



We discuss the galaxy-property functions in detail in lLauer et al. 



them briefly here. For the galaxy luminosity function we use the iBlanton et al.l (120031 ) Sloan 



(2007), but summarize 



Digital Sky Survey (SDSS) lu minosity function, t ransfo rmed to the l/-band and redshift 



z = 0.1, augmented with the iPostman fc Lauerl (119951 ) luminosity function of brightest 



cluster galaxies, which appear to be undercounted in the SDSS. We note that the Blanton et 
al. function refers to total rather than bulge luminosity in SO and spiral galaxies; however, 
this difference is less important at the bright end of the luminosity function, where the 
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selection effects are strongest. For the velocity-dispersion functio n we use the ISheth et al. 



( 120031 ) SDSS results, augmented at the high-cx end as prescribed by lBernardi et al. (Eoojto 
account for an artificial high-a cut-off in Sheth et al. Both functions are shown in Figure [TJ 

The function f(s) should be the least-squares fit of log 10 M. to log 10 L or log 10 cr in a 
sam ple selected by galaxy p roperties. The M. — a relation f(s) is based on the galaxy sample 



sec 



from Tremaine et al. 



Lauer et al. 



q> 

)2) 



2002) augmented by a few galaxi es with more recent M. determinations 



20071 ) . In contrast to the treatment in iTremaine eraD (booi ). which treated 



M, and a symmetrically, the appropriate treatment for our purposes is a least-squares fit of 
log 10 M. on log 10 a, which gives: 



log lo (M./M ) = (4.13 ± 0.32) log 10 (a/200 km s" 1 ) + 8.29 ± 0.07, 



(4) 



for H = 70 km s -1 Mpc -1 (which we will use throughout this pa per). We discuss the 
M. — L relation in detail in lLauer et al.l (120071 ); in brief we use the lHaring &: Rixl (120041 ) 
relation between M. and galaxy mass, transformed back to luminosity by adopting the 
m ass-to-light ratio M[Ly — 6 x io^ ao92 ( M ^+ 22 ) M / X , based on the M/L estimates given 



in 



Gebhardt et al.l (120031 ). A least-squares fit of M. on My for galaxies with My < —19 



gives: 



log lo (M./M ) = (1.32 ±0.14)(-M v -22)/2.5 + 8.67 ±0.09. 



(5) 



Equations (Jlj) and §5§ have the form 



f(s) = a + bs, 



(6) 



where a and b are constants. lLauer et al.l (120071 ) show that the M. — a and M. — L relations 
make different predictions for M. in luminous or high-dispersion galaxies, which have lower 
a values for a given luminosity than is implied by eliminating M, from equations (jlj) and 
flSJ). Thus at least one of the M. — a or M. — L relations must curve away from the linear 
fit ([6]) in the most luminous galaxies. We are not concerned with this issue here, and simply 
show the different selection biases that can result from using a or L as the galaxy property 
that correlates with M., assuming that the linear relation (jH]) holds at all masses. However, 
the reader should bear in mind that neither the M. — L or M. — a relation is well determined 
at high galaxy masses, where the potential selection bias is most important. 

The scatter function, h [/i — f(s)] , is poorly known at best. For this analysis, we assume 
that at any galaxy L or a, log 10 M. is described by a normal distribution about f(s), with 
cosmic scatter a M . Thus 



h [a* - f( s )} 4* 



7T(7„ 



exp 



\jj - /(a) 



(7) 
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Again, we emphasize that <r„ does not embody any observational errors i n the determ i nation 



of M.; it represents the intrinsic spread in M. at any galaxy property. iNovak et al.l (120061 ) 
conclude that only an upper limit to is presently known for either the M. — a or M. — L 
relation, both because of the small sample of reliable M. determinations, and because of 
uncertainties in the observational errors in M., which must be accurately measured to isolate 
the contribution of to the total residuals. For the present analysis we will explore the 
sensitivity of the selection biases to cr M on the assumption that < 0.3 for dispersion a and 
cr M < 0.5 for luminosity L. 

Apart from the poor knowledge of <t m , there are at least two other major uncertainties in 
the function h: (i) a more general treatment would allow for the possibility that it varies with 
s; (ii) there is little or no justification of the assumed normal form. As will be shown below, 
for the most massive black holes or galaxies, the selection bias may depend on the form that 
h [fi — f(s)} takes at several standard deviations away from the mean, yet an observational 
determination of the form of h in this region would require a sample of well-determined 
black-hole masses several orders of magnitude larger than is presently available. Thus at 
the highest black-hole masses, the selection bias is likely to depend sensitively on knowledge 
that is not presently at hand, so it is not possible to apply reliable corrections for this bias. 

With these caveats, in Figure [2] we show estimates of z/(/i, s) as modeled by equation 
(T5]) for each of the M. — a and M, — L relations, and two assumed values of a M . Projection 
of z/(/i, s) onto the /z-axis produces the mass function of black holes, 

$.(//) = j v{n,s)ds (8) 

h[fi-f(s)]g(s)ds. (9) 
The importance of cosmic scatter for converting the luminosi ty or dispersion function g(s) to 



2004; 


Tundo et al. 


2006; 


Lauer et al. 


2007) 



cosmic scatter implies that the most massive BHs are often hosted by modest galaxies that 
a priori would not be expected to harbor BHs of high mass. To illustrate this point, Figure 
[3] shows $,(/i) for four different versions of z/(/x, s). The salient feature is that as the cosmic 
scatter increases, the contribution to the density of the most massive black holes from the 
wings of the scatter function h overwhelms the "native" population of massive black holes 
harbored by the galaxies with the largest values of property s. A second illustration of this 
point occurs in the top two panels in Figure HJ which show the probability distribution of L 
and a for M. = 1O 1O M : note that the solid curves, corresponding to a M = 0.5 for L and 0.3 
for a, have a prominent shoulder to the right of the peak, the peak arising from low-luminosity 
or low-dispersion galaxies, and the shoulder from high luminosities and dispersions. 
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The source of the selection bias can be seen more directly in Figure [21 which shows 
z/(/i, s) for each of the M. — o and M. — L relations, and two assumed values of a M . Only for 
0^ = is there an exact relation between M. and o or L. Our model for z/(/x, s) — which may 
not be correct, but is consistent with the data — implies that the density contours in Figure 
[2] are symmetric about the mean relationship ridge-lines in the vertical direction; in other 
words, the conditional probability of ft given s is symmetric about \i = f(s). On the other 
hand, the conditional probability of s given /i, 

is not symmetric about s = f (ft), as seen in Figure HI which shows the distribution of s 
at selected values of M.. As M. and o^ increase, the distribution of s for a given M, moves 
further and further away from the value implied by the inverse of the M. relations, = 
(ft — a)/b, where a and b are the coefficients given in equation (jSJ). For M. > 1O 9 M , 
generally falls well out in the wings of the distribution of s for both L and o, particularly for 
the larger values of o^. The galaxies with s = f~ 1 (fi) for the highest mass black holes are so 
far down in the step cutoffs of the L and o distribution functions that they are completely 
overwhelmed by the population of galaxies of modest mass that harbor high mass black holes 
as statistical outliers. 

The mean value of s at a given M. is 

f fsv(fi,s)ds f sg(s)h[fi- f(s)] ds 

= / s u(s\fi) ds = ^— — - = f ■ (11) 

J Jv(fi,s)ds J g(s)h[fi- f(s)\ ds 

This is shown as the red lines in Figure [21 The mean of log 10 L or log 10 a at a given M. is 
offset from which is just the s location of the M. relation ridgelines shown in the 

Figure. In the presence of cosmic scatter, the mean s of a sample of galaxies that host a 
black hole of given ft is different from the mean fi hosted by a sample of galaxies of given s. 
This difference is the source of the selection bias. 



3. Illustration of Selection Biases 

Selection bias typically occurs because the galaxy samples used to probe the M. relations 
at cosmological distances are not selected by galaxy L or o, but by the visibility of their 
AGNs. For most of the discussion in this section we shall assume that AGN luminosity does 
not depend directly on L, o, or any other galaxy property unrelated to M,. This model is 
appropriate if, for example, the probability that the AGN associated with a black hole of 
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mass M. has luminosity Lagn is given by 

* = where f «*)te=l (12) 

and A = log 10 I/agn and /x = log 10 M.. The physical content of this assumption is that 
the luminosity history of an AGN is determined by the black-hole mass and scales with the 
Eddington luminosity, which is proportional to M.. 



3.1. Bias in a Luminosity-Limited Survey 

We first consider the bias that occurs in a sample in which (i) all of the objects are in a 
narrow redshift range; (ii) the survey contains all AGNs brighter than a given flux. In this 
case the probability that an AGN in the relevant redshift range is accepted in the survey 
depends only on its luminosity Lagn, and by equation (Tl2|) this in turn depends only on 
its black-hole mass fi. Thus the probability distribution of a galaxy property s at a given 
value of (jl is not biased by the selection, and is simply v(s\n) (eq. [TU]). Some examples are 
shown in Figure HJ which plots the probability distribution of s — for several values 

of black-hole mass M. and cosmic scatter a^. Here is the inverse of the M. — L or 

M. — a relation (eq. [2]). The mean values of the distributions (s)^ are displaced to values 
lower than in most of the examples shown. As expected, the width and offset of the 

distributions are larger for larger values of <r M . Also as expected, the offsets are generally 
larger for more massive black holes, (s)^ is shown as the red lines in Figure [2] the offset of 
these lines from the relation fi = f(s) demonstrates the selection bias. 

In this model, the selection bias can be corrected: an assumed form for the probability 
distribution z/(/i, s) such as fl2]) with redshift-dependent parameters can be used to compute 
v(s\n), which can then be fitted to the distribution of galaxy property s in the sample at a 
given value of black-hole mass \i (separate issues, discussed in §H are whether this model for 
the selection effects is realistic and whether the assumed form of u(fi, s) is correct). However, 
most papers in the literature have not taken this approach, so it is worthwhile to estimate 
the biases that might be introduced by using simpler statistics to estimate the evolution in 
the M. relations. 

Figure [5] shows the object selection bias As = — (s) M as a function of M, and 

a p. In this example we have expressed the bias in terms of galaxy properties, since we 
have selected by M.; however, it may also be represented as Alog 10 M. = bAs, where b is 
the slope given in equation ([6]). Note that the bias is not always a monotonic function of 
black-hole mass, nor does it always have the same sign. In this simple model, the bias at 
black-hole masses of 1O 9 M is Alog 10 M. = 0.4 in the M. — a relation for = 0.3, and 
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Alog 10 M. = 0.7 in the M. — L relation for = 0.5 — and even larger at larger black-hole 
masses. The bias is smaller at lower L or a, but we caution that the near-zero bias in the 
M. — L relation at sma ll M. is an artifact of our assumption that the luminosity function 



has the ISchechterl (119761 ) form at low luminosities (see below), and this may not be true for 



the early-type galaxy components that are believed to host the black holes. 

Figure [5] also shows that a significant portion of the bias comes from galaxies with 
\s — > 2a ^ for the larger values of a^. Indeed for M, k 5 x 10 9 M Q , nearly half of 

the bias comes from such galaxies for all but the smallest shown. Thus, not only is the 
bias sensitive to a M , it also depends on the shape of the wings of the error distribution. If, 
as is likely, the wings are more extended than in our assumed log-normal distribution, the 
bias will be even larger. 

Many of the features in these plots can be understood analytically. From the definition 
of (s)^ (equation [TTj) and equation (Q, we can write 

Jxh(x)g[b~\^-a-x)}dx 
+ {S) " f h(x)g\b-iQi -a-x)]dx' 1 ' 

If the cosmic scatter o 2 = J dxx 2 h(x) is not too large, we can evaluate this by expanding 
g(s) in a Taylor series, 

dlng(s 



a + &<*)„ = a* + *2 — +0(aj). (14) 

L as ls=(fi-a)/b 

This result shows that (i) the bias is proportional to a 2 if is not too large; (ii) the bias 
in (s)^ is positive if g(s) declines with s and proportional to dlng(s)/ds — this is why the 
bias is large and negative for the most luminous or high-dispersion galaxies; (iii) the bias is 
near zero if g(s) is constant, which corresponds to a luminosity function dn oc dL/L when 
s = log 10 L. The bias in the left panel of Figure [5] is small for low-luminosity galaxies because 
the assumed luminosity function is close to this form. 

A different statistical method, which more closely parallels the approach used in most 
observational papers, is simply to estimate the average value of the difference between \i and 
the prediction of the fi — s relation ([6]) for all the galaxies in the sample, 

A log 10 M. = {fi- f(s)) = { f i)-a- b(s), (15) 

Let us assume that the survey flux limit at the given redshift corresponds to a luminosity 
Lq with Ao = log 10 Lo- Then 

A log M - dX ^ g ^ ds I^~ a - hs )K^ ~ a ~ bs)jj(X - /j) d\i 
10 * dX J g(s) ds J h(n - a - bs)ip(X - fi) d/j, 
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where x = ji — a — bs, and 



dX J g(s) ds J xh(x)?p(X — x — a — bs) dx 
dX J g(s) ds J h(x)ip(X — x — a — bs) dx 
dX J xh(x)(f)(X — x) dx 
dX J h(x)(j)(X — x ) dx 



(17) 



Hv)= / g{s)^{y - a - bs) ds. (19) 



The luminosity function (number per unit volume) of AGNs is dn = *f?(X)dX, where 

\I/(A) = / g(s) ds / h([i — a — bs)tjj(X — fi) dfi = / h(x)(f)(X — x) dx, (20) 



which is just the inner integral of the denominator of ffTS]) . To evaluate the inner integral in 
the numerator of equation (fT8l) . note that the functional form of h(x) (eq. [7J implies that 
h'(x) = —xh(x)/a^, so 

x h(x)(f)(X — x) dx = — a" 2 . J h'{x)(j){X — x) dx. (21) 

Integrating the right side of the equation by parts, and noting that lim^oo (f>(y) = 0, gives 

x h(x)(f)(X — x) dx = a 2 h(x)—(p(X — x)dx (22) 

p ./ ox 



»dX 

Thus 



-u\ J h(x)-^-(f)(X-x)dx (23) 
-a 2 ^- [ h(x)(j)(X - x) dx. (24) 



Alog 10 M. = (25) 



Remarkably, the result depends only on the directly observable luminosity function ^(A) 
for the type of AGN targeted in the survey, and is independent of assumptions about the 
luminosity history %p(X — fi) or the distribution of host galaxy properties g{s). Also, in 
contrast to the analogous result (fT4|) . equation (1251) is not just the first term in a Taylor 
series in a 2 but valid for all values of a 2 ,, no matter how large. 

To illus t rate t he application of this result, we adopt the quasar luminosity function from 



Bovle et all (120031 ) . 

*( A ) = 10 -(l+a)(A-A*) + 10 _(i +/3 )(A-A*)> ( 26 ) 
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with a = —3.4 and (3 = —1.6, where A* is the "break" luminosity. The corresponding bias 
A log 10 M. is shown as a function of the lower luminosity limit Ao — A* in Figure [6j Note 
that the bias in Figure M becomes smaller but does not vanish as the survey goes to fainter 
and fainter luminosity limits: asymptotically, Alog 10 M. — ► —(1 + /3)lnl0 which equals 
0.12(a^/0.3) 2 . 



3.2. Bias in a Flux-Limited Survey 

As a second example, we examine the bias in a flux-limited sample of AGN. We assume 
that the survey galaxies are distributed uniformly in Euclidean space, that the probability 
distribution of AGN luminosities is given by equation (fT2|) . and that the survey contains all 
galaxies with flux / = Lagn/^ 2 exceeding some limiting flux f . 

The bias is then given by 

A log M — / ^ A / 9{s)nV{\) ds j{fi - a - bs)h(fi - a - bs)ip{X - /j) dji 
J dX J g(s)nV(X) ds f h(n — a — bs)ip(X — jj) djj, 

where ^(A) = |Afi(LAGN//o) 3 ^ 2 oc 10 3A//2 is the volume within which an AGN of luminosity 
A can be detected, AQ is the solid angle covered by the survey, and n is the number density 
of galaxies. We have 

/ dX jg{s)10 3X / 2 ds Jin- a- bs)h{y -a- fe)^>(A - fi) dfj, 
gl ° " / dX J 3( S )10 3A /2 ds j h{fx - a - bs)i>(X - fi) dfi 1 J 

/ 10 3A / 2 dX f g(s) ds fxh(x)iP(X — x — a — bs)dx 



f 10 3A /2 dX f g(s) f h(x)il>(X -x-a-bs)dx 
J 10 3A / 2 dXjx h(x)<f)(X - x) dx 
J 10 3A / 2 dX J h{x)<p{X -x)dx ' 



(30) 



where <fi(y) is defined by equation (1191) . To evaluate the inner integral in the numerator of 
equation fl30l we use equation (T2~4"l) . which yields 

, f 10 3A / 2 dX± [ h(xU(X -x)dx 
Alog 10 M. = vf \Z\ ' (31) 



/10 3A / 2 dX f h(x)4>(X - x) dx 



3 In 10 



< (32) 



where the last line follows from an integration of the numerator by parts. 
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3.3. Selection Effects That Depend Only on Galaxy Properties 

The results of the previous subsections depend on the assumption that the probability 
distribution of AGN luminosity is determined by the black-hole mass and not the galaxy 
properties (eq. [T2l . A foil to this is to assume that the luminosity is determined by the 
galaxy properties and not the black-hole mass. Thus equation (IP21) is replaced by 

dp= X (X-ks)dX where =1 (33) 

where k is a constant. In this case A log 10 M. = in any flux-limited or luminosity-limited 



sample (ISalviander et al.ll2006l ). 



4. Correcting for Selection Bias May be Extremely Difficult 

Given the example bias calculations shown in the previous section, it may be tempting 
to conclude that similar bias corrections may be estimated and applied after the fact to 
existing surveys to determine the M. — a or M. — L relations. We believe that in general such 
corrections are difficult or impossible to apply reliably to current surveys, for the following 
reasons. 



4.1. The Error Model is Poorly Known 



In the examples above, we assumed that at any galaxy property s, \i = log 10 M. followed 
a normal distribution about the mean M. — s relation /i = f(s), with standard deviation 
characterized by a cosmic scatter cr M . We saw that the bias is very sensitive to (in most 
examples oc cr^) yet this parameter is poorly known for either the M. — a or M. — L relations, 
in part because an a ccurate estimate of <7 M requires knowing the observational errors in the 
M. determinations. iNovak et al.l (120061 ) have studied this problem and conclude only that 



Qiu 5, 0-3 for the M , — a relation and oy, < 0.5 for the M. — L relation. Moreover, the 



Novak et al.l (120061 ) analysis assumes, as we do, that cr M is constant over the parameter range 
of the M. relations; the sample of galaxies with well-determined M. is simply too small to 
explore the possibility that it is not. 

Even if an accurate estimate of <7 M were available, the functional form of the distribution 
of \i — f(s) is unknown. The assumption of a normal distribution is an obvious first step, yet 
the steep falloff of the L or a distribution at large values means that the selection bias for 
large black-hole masses is sensitive to the wings of this distribution, where \x — f(s) > 2cr M . 
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The assumed normal distribution is likely to underestimate the selection bias since more 
realistic distributions have fatter tails. Given that the local black-hole sample is so small 
that is poorly determined, the sample of galaxies with good M. determinations would 
have to be several orders of magnitudes larger before an accurate form for the scatter function 
could be determined. 



4.2. The M. Relationships are Poorly Known 



The M, relations are imperfectly known at z = 0, and in particular the range of L 
or a that is used to determine them is rather limited. For example, i n the sample of 31 
local galaxies with measured black-hole masses used by iTremaine et al.l (120021 ) to determine 
the M. — a relation, the interquartile ranges in L and a are only factors of 5.3 and 1.6 
respectively. Thus, the calculation of the selection biases for M. > 10 9 M Q must contend 
with the fact tha t there are only a few black hole s observed in this mass range in the local 
sample (4 in the ITremaine et al.l 120021 sample) . lLauer et al.l (120071 ) show that the L — a 
relation must be curved (in logarithmic coordinates) in the sense that a appears to increase 
only slowly (if at all) with L for the most luminous galaxies. This implies that the present 
log-linear M. — a and M, — L relation cannot be consistent at high galaxy luminosity. The 
joint M. — s distributions shown in Figure [2] are thus based on extrapolated estimates at the 
highest M. values that are likely to change if more and better determinations of black-hole 
masses M. > 10 9 M Q become available. 



4.3. What Determines AGN Luminosity? 

The heart of our analysis of selection bias in § §3.11 and 13.21 is the assumption that the 
black hole determines the AGN luminosity, or more precisely that the galaxy properties do 
not. This is clearly true if, for example, (i) black holes radiate either at the Eddington 
luminosity LEdd or not at all, and (ii) the bolometric correction is independent of galaxy 
properties. This assumption is also true in more general circumstances, for example, (i) can 
be replaced by the weaker assumption that the probability that a black hole is radiating 
at some Lagn depends only on the ratio L^cm/L-fm (eg . [T21) . This assumption may be 



approximately correct: for example, iHopkins et al.l (120061 ) find that their simulations are 



well-fit by a model in which the probability distribution of AGN luminosities depends only on 
ZvAGN/^peak where Lpeak is the peak luminosity of the AGN and L pea k oc £ E dd- Nevertheless, if 
galaxy properties do affect the AGN luminosity — most likely by determining the feeding rate 
of matter onto the black hole at luminosities L <C L-^dd — then the biases will be different 
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from those presented in § §3.11 and 13. 21 In the limiting case that the AGN luminosity is 
determined entirely by galaxy properties rather than black-hole mass, as might be plausible 
for low-luminosity AGNs, there is no selection bias in a flux- or luminosity-limited survey 

4.4. Did the Galaxy Make the Black Hole, or the Black Hole the Galaxy? 

We have assumed that the joint probability distribution of black-hole mass \i and galaxy 
property s can be modeled by equation (j2J). This corresponds to a physical model in which 
the galaxy property is the independent variable, and the black-hole mass is determined by 
the properties of the galaxy through the cosmic-scatter function h(y) and the M. relation 
f(s). This assumption is convenient for modeling the local sample of galaxies with M, 
determinations, which was selected by galaxy properties, plus the ready knowledge of g(s), 
the volume distribution of s. The alternative model provided by equation ([3]), in which the 
black-hole mass is the independent variable and the galaxy property is determined by the 
black hole through the function p(y), is much more difficult to fit to the local observations as 
we have little direct knowledge of the black- hole mass function. However, our present 

understanding of black-hole and galaxy formation does not allow us to say which model, if 
either, is correct. 

4.5. Survey Selection Effects May be Poorly Known 

Most of our estimates of selection bias have been based on the assumption that surveys 
are complete to a given AGN luminosity or flux, but this is an oversimplification. Any survey 
that is based on measurements of the luminosity of the host galaxy requires that the host is 
bright enough and large enough to be separated from the AGN. Any survey that uses the 
velocity dispersion measured from stellar absorption lines requires that the absorption lines 
are strong enough, and that AGN emission lines in the vicinity of the absorption lines are 
weak enough, to allow a reliable dispersion measurement. At the opposite extreme, surveys 
of low-luminosity AGNs require that the emission lines are strong enough to be detected 
against the continuum flux from the galaxy. In such studies the selection depends in a 
complex way on properties such as the ratio of AGN to galaxy luminosity, and on galaxy 
properties other than L or a, such as the effective radius or central surface brightness. 
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5. Bias in Existing Surveys 



We have intended this paper mainly as a planning guide for future surveys to explore 
the evolution of the M. relation, rather than as a critique of existing surveys. The following 
brief comments on existing surveys are mostly intended to illustrate the application and 
impact of the estimates of selection bias that we have made in earlier sections. 

Several authors have determined the M, — a relation for nearby AGNs, and compared 
this to the local relation for inactive galaxies. Using seven Seyfert 1 g alaxies with black- 
hole masses measured by reverberation mapping, iGebhardt et all (12000br) find Alog i n M. = 
—0.21 ± 0.13. For 16 AGNs with reverberation-based black-hole masses lOnken et all (120041 ) 
find Alog 10 M. = —0.26 ± 0.15 (although Onken et al. interpret their result in terms of 
calibrating the reverberatio n-mapping method rather than an offset in the M. — a relation). 
Barth. Greene. &: Hoi (120051 ) have compared the black-hole masses in dwarf Seyfert 1 galaxies 
to an extrapolation of the M. — a relation to lower dispersions; using black-hole masses 
based on the photoionization method calibrated by the standard "isotropic" formula, they 
find A log 10 M. = —0.04 (although t hey interpret t heir r esults with Onken et al.'s calibration 
and thus find A log 10 M. = +0.23). iGreene fc Hoi (120061 ) have measured the M. — a relation 
in 88 nearby AG N using the photoion ization method to determine M. — the sample is much 
larger than that of lOnken et all (120041 ) but the photoionization mass estimates are less direct 
than those from reverberation mapping. They find Alog 10 M. = —0.21 ± 0.06. 



Similarly, lLabita et al.l (120061 ) have determined the M. — L relation from a sample of 29 
quasars with z < 0.6; they determine black-hole masses using the photoionization method 
and luminosities using HST photometry. They find Alog 10 M. = —0.36, although they 
interpret their result in terms of recalibrating the photoionization method. 



A com mon feature of t he IGebhardt et al.l (l2000bl ) , lOnken et al.l (120041 ) , iBarth. Greene, fc Ho 



( 120051 ). and I Greene fc Hoi (120061 ) samples is that their black-hole masses are relatively mod- 
est: generally M. < 1O 8 M and often more than an order of magnitude smaller. Selection 
bias occurs at all M., depending on the local slopes of the L and a distributions, but these 
samples are safely removed from the very strong biases seen at large values of L, a, and M., 
where the distribution functions are falling rapidly. Although we have emphasized selection 
bias in the context of low- and high-redshift samples with similar black-hole masses, one 
selected by galaxy properties and one by black-hole mass, selection bias can also arise if 
the two samples are selected in the same way but are centered on different black-hole mass 
ranges. 



We furt her note that sel e ction e ffects in these samples are difficult to model — for ex- 



ample, 



U as 



Gebhardt et al.1 (l2000bl ) point out, the AGN variability timescale depends on 
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luminosity, and if the timescale is too long, reverberation mapping is impractical; (ii) to be 
able to measure the velocity dispersion requires some minimum ratio of the luminosity of 
the host galaxy to the luminosity of the AGN. If we make the simplest possible assumption, 
that the samples of nearby galaxies are complete but flux-limited, then equation (1521) implies 
Alog 10 M. = 0.31(V0.3) 2 . 

Woo et all 6x3), building on the initial work of iTreu et al.l (120041 ). examine a sample 
of 14 Seyfert 1 galaxies at z = 0.36 ± 0.01, measuring the velocity dispersions from stellar 
absorption lines and the black-hole masses by the photoionization method. Compared to the 
local sample of inactive galaxies with measured black-hole masses, they find A log 10 M. = 
0.62 ± 0.10 ± 0.25; note that this result is based on a calibration of the photoionization 
masses that assumes that the lOnken et al.l (120041 ) sample should have exactly the same 
M. — a relation as inac tive galaxies, which need not be the case if selection bias is accounted 
for. iWoo et al.l (120061 ) argue that their offset A log 10 M. is robust because the black-hole 
masses are measured by the same technique, with the same calibration, in the local and high- 
redshift samples, bu t the selection bias in the two samples is likel y to be quit e differ ent; the 
Onken et al.l (120041 ) sample is more nearly flux-limited, while the IWoo et al.l (120061 ) sample 
is more nearly luminosity-limited, and the selection bias in these two cases is quite different 
(eqs. [32]and[25]). A useful next step in assessing t he bias would be to determine the luminosity 
function for AGN selected by the criteria used by lWoo et al.l (120061 ). for use in equation (125]) . 



It may also be significant that the typical M, Ri4x 1O 7 M in the lOnken et al.l (120041 ) 



sam ple is roughly an order of magnitude less massive than the typical M. « 3 x 1O 8 M in 



the IWoo et al.l (120061 ) sample. As is evident in Figure [5j the bias is a strong function of M., 
and indeed the bias in the M. — a relation changes sign at M , ~ 10 8 Thus we would 



expect that selection bias might affect the r esults even if the lOnken et al.l (120041 ) sample 
were chosen in precisely the same way as the IWoo et al.l (120061 ) sample. 



Figure [5] predicts a bias of Alo g in M. ~ 2 at M . « 3 x 1O 8 M for a M = 0.30. This 
is only 40% of the offset claimed by IWoo et al.l (I2006 ). Th us, it appears that selection bias 
cannot explain all or even most of the IWoo et al.l (j2006l ) Alog 10 M.; however, given the 
uncertainties discussed in the previous section, the possibility that the correction may be as 
large as the claimed offset cannot be ruled out. 



Peng et al.l (J2006J) find A log 10 M. ~ —0.1 compared to the local M, — L relation, using a 
sample of 11 quasars at z ~ 2; the black- hole masses are determined with the photoionization 
method. If we assume that the velocity dispersion of the host galaxy does not change between 
z ~ 2 and z = 0, and that the luminosity of the host fades by 1-2 magnitudes as predicted by 
passive stellar evolution, then this result corresponds to Alog 10 M. ~ 0.4-0.8 as measured 
by the M, — a relation. However, the potential selection bias is rather large for this sample, 
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partly because the cosmic scatter in the M. — L relation may be larger than in the M. — a 
relation. If the sample is luminosity limited, the model shown in the left panel of Figure 
predicts Alog 10 M. = (0.3-0.7) for = 0.5 for the range of black-hole masses in the 
sample, and the flux-limited model predicts Alog 10 M. = 0.58 for = 0.5 (note that our 
assumption of a Euclidean metric is not valid at z ~ 2, and should be generalized to a realistic 
cosmological model). A further complication is that their method requires measuring the 
bulge luminosity in the presence of the bright AGN, so imposes a selection on the ratio 
of AGN to galax y luminosity. Without an accurate estimate of the selection bias for the 



Peng et all (120061 ) sample, their measurement of A log 10 M. does not provide clear evidence 



for rapid evolution in the ratio of black-hole mass to stellar mass in the galaxy. 



Salviander et all (120061 ) find that galaxies of a given dispersion at z ~ 1 have black-hole 
masses that are larger by A log 10 M. ~ 0.45 than at z ~ 0. In contrast to most other studies, 
they carefully consider the effects of selection bias using Monte-Carlo models and conclude 
that for 0^ = 0.3 selection bias increases log 10 M. by ~ 0.1 at low redshift and 0.2 at z ~ 1, 
contributing a net offset of 0.1. They also argue that scatter in the relationship between the 
width of the [O III] emission line and the true stellar dispersion contributes an additional 
offset of ~ 0.15, so after correcting for these biases Alog 10 M. ~ 0.2. Their estimates of 
selection bias are similar to, but somewhat smaller than, the bias that we estimate for their 
sample from equations (|25j) and (|32"1) . assuming that the low-redshift sample is flux-limited 
and the high-redshift sample is luminosity-limited; we believe our estimates are more robust 
because they do not depend on an assumed form for the distribution ^(A — /i) of AGN 
luminosities at a given black-hole mass. 



Finally, we note that theoretical models of the evolution of the M.— a relation ([Robertson et al 



20061 ) predict Alog 10 M. = —0.18 at z = 2. Strong selection bias could reduce or eliminate 
the difference between this prediction and the positive values of A log 10 M, found at high 
redshift in most of the above studies. 



6. Recapitulation 

We have described a selection bias that affects investigations into the cosmological evo- 
lution of the M. — o or M. — L relations. We summarize the main points as follows: 

• Cosmic scatter in the M, — o or M, — L relations means that a galaxy of given L or 
a will host black holes with a range of M. values. 

• The steep decline in the luminosity function for the most luminous galaxies means that 
the rare occurrence of high-mass black holes in numerous "modest" galaxies overwhelms 
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the frequent occurrence of black holes of similar mass in the rare galaxies with very 
high luminosity. A similar conclusion applies to dispersion instead of luminosity. 

Selection of a sample of black holes hosted by inactive galaxies at low redshift explores 
the distribution of black-hole masses M. that a galaxy will host, given its luminosity 
L or dispersion a. Selection of a sample by AGN luminosity at high redshift explores 
the distribution of L or o at a given M. (if, as is likely, the distribution of AGN lumi- 
nosity is determined mainly by the black- hole mass rather than the galaxy properties). 
Confusing the two distributions may create a false signature of evolution. 

The selection bias is substantial for the current estimates of the cosmic scatter in the 
M, — L or M. — a relations, particularly for galaxies containing the most massive 
black holes. Neither the root-mean-square cosmic scatter nor the shape of the scatter 
function is known accurately and hence the bias cannot be reliably corrected for, given 
our present and even foreseeable state of knowledge. 

Our rough estimates of selection bias as applied to a number of evolutionary studies 
show that the potential bias may be larger t han many of the A log -, n M, or As val 



ues ob served. With the notable exceptions of iFine et al.l (120061 ) and ISalviander et al. 



(120061 ). it appears that most studies of the evolution of the M, relations have not con- 
sidered the effects of selection bias. In most cases, accounting properly for selection 
bias is likely to reduce the the evolution in the M, relations that has been observed in 
recent surveys. 

The only way to avoid selection bias is to choose high- and low-redshift galaxy sam- 
ples using precisely defined, objective criteria that are precisely the same for the two 
samples. To our knowledge, this has not yet been done. 



The work was initiated during a Lorentz Center workshop attended by the authors. We 
thank the staff of the Lorentz Center and Dr. Tim de Zeeuw for hosting us. We thank Drs. 
Karl Gebhardt, Richard Green, and Nadia Zakamska for useful conversations. 
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Fig. 1. — The number density of galaxies as a function of luminosity (left) and velocity dis- 
persion (right), assuming H = 70 km s -1 Mpc -1 . The luminosity function is based on data 
fro m the Sloan Digital Sky S urvey (IBlanton et al.ll2003l ). augmented at large luminosities by 
the iPostman fc Lauerl (119951 ) sa mple of brightest c luster galaxies. The velocity-dispersion 
function is based on SPSS data fjSheth et al.l 120031). augmented at large dispersions by the 
sample identified by Bernardi et ali ( 20od ). See Lauer et al. ( 2007 ) for details. 
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Fig. 2. — The joint probability distribution of galaxy luminosity L (top panels) or stellar 
velocity dispersion a (bottom panels) and black- hole mass, M.. The solid, straight, black 
line gives the mean M. — L or M. — a relation (eqs. H]and[5]). The joint distribution is 
calculated by assuming that the black-hole mass has a log-normal distribution about the 
M. — L or M, — a relation, with the normalization provided by the galaxy luminosity or 
velocity-dispersion functions shown in Figure [TJ The adopted dispersions in log 10 M. are 
ay, = 0.25, 0.50 for the M, — L relation or = 0.15, 0.30 for the M. — a relation. The 
contours are arbitrary, but are spaced in increments of 0.5 dex in density. The curved red 
lines give the mean My or log 10 (<j) as a function of M.. The displacement of this line from 
the mean relations gives the bias at a given M.. 
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Fig. 3. — The effects of cosmic scatter in the M. — a or M. — L relations on the number 
density of black holes as a function of M, . The left panel is based on the M. — L relation of 
equation (jSJ) combined with the luminosity function of Figure [TJ The right panel is based on 
the M. — a relation of equation (j3J) combined with the velocity-dispersion function of Figure 
[U Solid lines show the results assuming no scatter in the M. relations, while the dashed 
lines show the effects of increasing amounts of cosmic scatter cx M . The dotted and dashed 
lines show the effects of = 0.25, 0.50 for the M. — L relationship and <r M = 0.15, 0.30 for 
the M. — u relationship. Cosmic scatter greatly enhances the volume density of black holes 
at the high-mass end. 
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Fig. 4. — The distribution of absolute magnitude My = —2.5 log 10 L + const or velocity 
dispersion log 10 cr at selected values of M. (log 10 M m /M Q is given in the upper right of each 
panel). These probability distributions are obtained by taking horizontal cuts through the 
joint distributions shown in Figure and are appropriate when the sample selection is on 
black-hole mass. Two distributions are shown with dashed and solid lines, corresponding 
respectively to ct m = 0.25, 0.50 for My on the left or = 0.15, 0.30 for log 10 a on the right. 
The horizontal coordinate is the difference between either My or log 10 a and the nominal 
values / _1 (s) determined from the M. — L or M. — a relations (j3J) and (jSJ). 
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Fig. 5. — The ratio of the nominal values determined from the M, — L or M, — a 

relations to (s) M , the mean of log 10 L or log 10 a at a given black-hole mass. This Figure 
quantifies the bias in a sample selected by M, as opposed to s. The blue curves are drawn 
for cr^ = 0.1, 0.2, . . . , 0.5 for log 10 L (left panel), and <r M = 0.1, 0.2, 0.3 for log 10 a (right panel). 
The cyan curves show the biases that result if the assumed normal distributions in log 10 M, 
(eq. Ej) are truncated at ±2a M . At high M. the bias depends not only on the amplitude of 
the cosmic scatter, but also on the form of the wings of the scatter function. 
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Fig. 6. — The bias in log 10 M. from the M. — a relation for inactive galaxies in a luminosity- 
limited survey of AGN. The horizontal axis is A — A* = log 10 (L /L*) where L is the limiting 
luminosity at the survey redshift and L* is the break luminosity of the AGN luminosity 
function (1261) . 



